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Abstract. We mainly investigate a coupled system of the generalized nonlinear Schrodinger 
equation and the Maxwell-Bloch equations which describes the wave propagation in an erbium- 
doped nonlinear fiber with higher-order effects including the forth-order dispersion and quin- 
tic non-Kerr nonlinearity. We derive the one-fold Darbox transformation of this system and 
construct the determinant representation of the n-fold Darboux transformation. Then the de¬ 
terminant representation of the nth new solutions p^'^\ which were generated from 

the known seed solutions (E, p, rf) is established through the n-fold Darboux transformation. 
The solutions p^'^\ provide the bright and dark breather solutions of this system. 

Furthermore, we construct the determinant representation of the nth-order bright and dark 
rogue waves by Taylor expansions and also discuss the hybrid solutions which are the nonlinear 
superposition of the rogue wave and breather solutions. 


PACS numbers: 42.65.Tg, 42.65.Sf, 05.45.Yv, 02.30.Ik. 

Keywords: Generalized nonlinear Schrodinger and Maxwell-Bloch system; Darboux transfor¬ 
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1. Introduction 

Recently, the long haul optical communication has attracted considerable interest of sci¬ 
entists around the world. But the efficiency of propagation of optical communication is still 
not very well. There are two important reasons, one is because of the dispersion, the other 
is due to the attenuation. In telecommunications, what we are interested is the variation of 
group velocity with frequency, because the absolute wave phase is often not important while 
the propagation of pulses is important. The dispersion makes the spread of the optical pulse 
temporally and may lead to the falling of the energy on the next bit slot. The dispersion is 
the linear effect for the propagation of optical pulse in fibers. The attenuation results from the 
optical losses which are the inherent feature of the optical fiber. The optical losses also cause 
the vanishing of the optical pulse due to the absorption and scattering . 

*Gorresonding Author. 
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What is important for the propagation of optical pulses in optical fibers is the nonlinear 
effect. The optical fiber behaves nonlinearity when the intensity of the optical pulse exceeds 
a certain threshold value. The most crucial effect is the self-phase modulation (SPM). While 
traveling in fibers, an optical pulse will induce a varying refractive index of the fiber due to 
the Kerr effect. This variation in refractive index will produce a phase shift in the pulse 
which leads to a change of the pulse’s frequency spectrum. The spectral broadening process 
of SPM can balance with the temporal compression due to the anomalous dispersion and 
reach an equilibrium state when the pulse is of adequate intensity. The resulting pulse is 
called an optical soliton [^. The possibility of the propagation of optical solitons which are 
governed by the nonlinear Schrodinger (NLS) equation was firstly introduced by Hasegawa 
and Tappert in 1973 [^|^. Mollenauer et al. observed experimental solitons in low-loss fiber 
in 1980 [^. Another momentous nonlinear effect is the self-induced transparency (SIT). The 
SIT means that the coherent absorption and re-emission of pulses make the two-level medium 
optically transparent to this wavelength when the energy difference between the two levels of the 
medium matches with the optical wavelength. The optical pulse can be amplified and reshaped 
by passing through an active zone doped with resonant atoms like the erbium. The resonant 
interaction can neutralize the optical losses in fibers. McCall and Hahn put forward SIT solitons 
in a two-level resonant system in 1967, which is usually described by the Maxwell-Bloch (MB) 
equations [^. 

Considering the optical fiber doped with resonant materials such as the erbium, which is 
governed by a coupled system of the NLS equation and the MB equations, the optical pulses 
satisfy both the NLS equation and the MB equations. It is necessary to amplify the optical 
pulses because of the decaying in the process of the propagation in fibers. In 1983 Maimistov 
and Manykin firstly proposed the NLS-MB system [^, after that Nakazawa and his cooperators 
observed SIT solitons in an erbium-doped silica fiber in 1991 and more studies about the 
NLS-MB can be seen in Refs. [9-12 . 


Mitschke and Mollenauer found that the observed solitons in experiments did not match 
with the theoretical properties [^ . The propagation equation approximated to a second-order 
dispersion (group velocity dispersion) with a cubic Kerr nonlinearity (SPM). This leads to 
the difference between the experimentally observed soliton and NLS soliton in normal fibers. 
The difference resulted from the additional perturbation of higher order effects such as the 
higher-order dispersion, self steepening and higher-order nonlinear effects. Therefore the same 
difference will appear, if one experimentally studies the NLS-MB soliton in erbium-doped optical 
fibers. So it is necessary to investigate the propagation equation with higher-order effects. The 
generalized nonlinear Schrodinger (GNLS) equation is derived by Porsezian and his partners 


14 . The coupled generalized nonlinear Schrodinger and Maxwell-Bloch(GNLS-MB) system is 
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described as 



Ez — i{Ett + 2\E\ E) + ir^Etttt + 8\E\ Ett + 2E‘^E^ 
+ 6E*E^ + 4\EtfE + GlEl'^E) + 2p, 

Pt = 2iujp + 2Ep, 

Pt = -{Ep* +pE*), 


tt 


(l.la) 

(1.1b) 

(1.1c) 


where subscripts z, t denote as partial derivatives with respect to the distance and time, the 
asterisk symbol as the complex conjugate, E as the normalized slowly varying amplitude of the 
complex held envelope, p = as the polarization, and p = \vi'^ — \v 2 f‘ as the population 
inversion with Vi and V 2 representing the wave functions of the two energy levels of the resonant 
atoms, u) as the frequency. The more crucial reason to study the above system in this paper is 
because of the existence of the quintic non-Kerr nonlinear term which is more signihcant than 
the cubic Kerr nonlinearity because the non-Kerr nonlinearity is responsible for the stability of 


localized solutions 17,18 


In recent years, comparing with solitons, the study on rogue waves in optics has also at¬ 
tracted considerable research due to their potential applications in different branches of physics. 
The study started from the pioneering measurement of Solli et ah by analyzing super-continuum 


generations in optical hbers 19 . The rogue waves appear from nowhere and disappear withont 


a trace, which is the description of the characteristics of rogue waves (^. The rogue wave 


occurs for the modulation instability (MI) |^- 25 . One of the possible generating mechanisms 
for rogue waves is the creation of breathers, then the larger rogue waves can be generated when 


two or more breathers collide 26 . The research on rogue waves has made many achievements 


among which Akhmediev has reported the recent progress in investigating optical rogue waves 
in Ref. |^ . 

To the best of our knowledge, there are few people to stndy the GNLS-MB system in the 
Eq. E3 so far. The solitons and breather solutions of the GNLS-MB system have been partly 


established in Ref. 16 , but the rogue waves of this system is still not reported by anybody. 


We will construct the determinant representation of the n-fold Darboux transformation of the 


GNLS-MB system, which is similar to the NLS-MB system 12 and H-MB system I28l-l30l . Then 


the nth-order rogue waves of the three optical fields will be given by determinants. Moreover, 
the p and p are found to be dark rogue waves. 

The paper is organized as follows. In Section the Lax pair of the GNLS-MB system 
is recalled, and we derive the one-fold Darboux transformation of the GNLS-MB system. In 
Section the determinant representation of the n-fold Darbonx transformation and formnlas of 
(£;N^ pN, ^N) are expressed. In Section]^ the bright and dark breather solutions are generated 
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from periodic seed solutions. In Section we construct the determinant representation of the 
nth-order rogue wave by applying Taylor expansions and discuss the effects of parameter r on 
the rogue wave solutions. Additionally, we discuss the hybrid solutions which are the nonlinear 
superposition of the rogue wave and breather solutions. Finally, we summarize the results in 
Section |6l 


2. Lax pair and the one-fold Darboux transformation 


In this section, we will derive the one-fold Darboux transformation of the GNLS-MB system 
in the Eq. El of which the Lax pair is 


Tt = [/T, (2.1a) 

= I/T, (2.1b) 


where 


with 


U 

V 


/^i(A; t, z)\ 
yT 2 (A; t, z)) ' 

—iXa^ + Uq, 

AV 4 + A" Eg + AV 2 + AEi + 1/0 + 

A -h C(j 


1 0 


= 


^0 -1 
1/4 = Sfrcrg, 


Uo = 


V 3 = - 8 tUo, 



Vo = 


Ei = 


En = 


—2i — 4fr \Ef —AirEt 
—AirE^ 2i-|-4ir lEl^ 

2TE*Et - 2 tEE; 

-2E* - At\E? E* -2TEt 


2E + At\E\^ E + 2TEtt 
-2TE*Et + 2 tEE* 


ail (I 12 
, 0-21 — 0 - 11 , 





-P 


(2.2a) 

(2.2b) 

(2.2c) 


Uii = i \Ef + Sir \Ef — ir \Et^ + iTE*Ett + irEE^^, 
ai2 = iEt + 6iT \E\^ Et + irEut, 

021 = + dir \Ef E^ + 

Here A is a complex parameter. 
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Basing on the Darboux transformation for the Ablowitz-Kaup-Newell-Segur(AKNS) sys¬ 
tem [^, we consider the following transformation of the Eq. (0. 


■irW = r,» = + s)», 


(2.3) 


where 


such that 


Ti = XI + S, 



^[ 1 ]^ = 
^[ 1 ]^ = 


(2.4a) 

(2.4b) 


Here depend on ew , , 7 ]^^^, A and they have the same form as U and V by replacing 

E, p, 7] by A'W, pX\ In order to make the Eq. ( |2.4[ ) invariant under the transformation 
(2.3), the Ti must satisfy 


ru + Tit/ = t/Wri, (2.5a) 

Ti, + TiV = V^^^Ti. (2.5b) 


By comparing the coefficient of A* {i 


0, 1, 2) on both sides of the Eq. (2.5a), we have 


EW =E + 2iSi2, S 21 = -S* 2 , (2.6a) 

Ai = [[/o, A]+z[(T 3, A]A. (2.6b) 


On the other hand, by multiplying by (A -|- cj) and comparing the coefficient of A* (i = 
0, 1, ..., 6) on both sides of the Eq. (2.5b|), we get 


eW=E + 2zAi 2, (2.7a) 

= (A - cjI)V-i(S - ul)-^ 

= TiU=_,E_irri|A=_,. (2.7b) 


Additionally, there are some constraints for elements of eI\^, for instance, (E )22 — 0, 

p\X\ = (e[\^)ii g R and (eI \^)^2 = (^ii) 2 i- The key step is to find the specific form of S expressed 
by the column solution of the Eq. ( |2.1[ ). Let 

S=-HKH-^, ( 2 . 8 ) 
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where 



f'i/iiXi; t, z) 'i/i{X 2 ; t, z)\ 
y^^ 2 (Ai; t, z) ^^ 2 (^ 2 ; t, z)J 


(2.9a) 


(2.9b) 


and det{H) 7 ^ 0, Ai and A 2 are complex constants. 

In order to satisfy the constraints of S and as mentioned above, we take the following 
constraints 


A2 = At, 

'^i{Xi,t,z) -^*{Xi, t, z)' 

d' 2 (Ai; t, z) d^t(Ai; t, z) 


H = 


( 2 . 10 a) 

( 2 . 10 b) 


So after taking the Eqs. (2.8) and (2.10) back into the Eq. (2.7), it results in the following 


one-fold Darboux transformation of the GNLS-MB system 


EW = E + 2z5i2, 

1 

det{Ti) 

1 


p\X\ = 


-2r^(Ti)n(ri)i2+P*(ri)i2(Ti)i2-p(Ti)n(ri)n]|A=-.., 


( 2 . 11 a) 

( 2 . 11 b) 


r^W = 


det{T\) 


[r?((Ti)n(ri)22 + (Ti)i2(Ti)2i) -p=^(ri)i2(Ti)22 +p(Ti)n(Ti)2i]|A=-u,. (2.11c) 


3. N-FOLD Darboux transformation for GNLS-MB system 


In this section, we will establish the determinant representation of the n-fold Darboux 


transformation for the GNLS-MB system as in Ref. 31 . For this purpose, we need to introduce 


2n eigenfunctions by = fk{Xk) = {fll) associated with an eigenvalue A^, and A^ 7 ^ A^ if 
k ^ m, where k = 1, 2, , 2n. Additionally, the eigenfunctions for distinct eigenvalues are 

linearly independent. 

According to the form of Ti, the n-fold Darboux transformation should have the form 
Tn = Tn{X) = X^I + + t 2 X^~^ + • • • + tn-\X -\- G, where tj{j = 1, 2 , ..., n) are 2 x 2 

matrices. From the following identity 


Tn{X', Ai, A 2 , A 3 , A 4 , . . . , A 2 n- 1 , X2n)\\=\kfk — 0 , (/c — 1 , 2 , ..., 2 n), (3.1) 

we can get the coefficients tj{j = 1 , 2 , ..., n) by the Cramer’s rule. Thus we obtain the 
determinant representation of the T„ in the following theorem. 


6 








Theorem 1. The n-fold Darboux transformation of the GNLS-MB system is Tn = T„(A) = 
A”'/ + ti\^~^ + + • • • + tn-i\ + tn, whose determinant representation is 


where 


W2n = 


(Tn)ll = 


(Tn)l2 = 


(Tn)21 = 


T -- 

- Tn{X; Ai, A2, 

A3, A4, . 

• ) A2n-1, 


f (Tn)ll 

V(^n)21 

hi 

/12 

Al/ll 

Al/12 

A?/ii 

a;/i2 ... 

ArVii 

/2I 

f 22 

A2/2I 

A2/22 

A2/21 

A2/22 

ArV2i 

/3I 

/32 

A3/3I 

A3/32 

Ai/31 

Ai/32 

ArV 3 i 

\yf2nl 

/ 2 n 2 

A 2 n/ 2 nl 

A 2 n/ 2 n 2 

1 — 1 

^ 2 nhn 2 ■ ■ 

A^-V 2 nl 

1 

0 

A 

0 

A 2 

0 

An -1 

fll 

/12 

Al/ll 

Al/12 

A?/n 

A?/i 2 . . . 

ArVii 

/2I 

f 22 

A2/21 

A2/22 

Ai/21 

A2/22 

ArV 2 i 

/3I 

/32 

A3/31 

A3/32 

Ai/31 

Ai/32 

ArV 3 i 

f 2 nl 

/ 2 n 2 

A 2 n/ 2 nl 

^ 2 nf 2 n 2 

^ 2 nhnl 

^ 2 nhn 2 ■ ■ ■ 

A^-V 2 nl 

0 

1 

0 

A 

0 

A 2 ... 

0 

fll 

/l 2 

Al/ll 

Al/12 

A?/ii 

A?/i 2 . . . 

ArVii 

/2I 

f 22 

A2/21 

A2/22 

A2/21 

A2/22 

ArV 2 i 

/3I 

/32 

A3/31 

A3/32 

Ai/31 

Ai/32 

ArV 3 i 

f 2 nl 

/ 2 n 2 

^ 2 nf 2 nl 

A 2 n/ 2 n 2 

1—1 

^ 2 nhn 2 ■ ■ ■ 

AL“V 2 nl 

1 

0 

A 

0 

A 2 

0 

An -1 

/ll 

/12 

Al/ll 

Al/12 

A?/ii 

Ai/12 

ArVii 

/2I 

/22 

A2/21 

A2/22 

Ai/21 

Ai/22 ... 

ArV 2 i 

/3I 

/32 

A3/31 

A3/32 

Ai/31 

Ai/32 

ArV 3 i 

f 2 nl 

/ 2 n 2 

^ 2 nf 2 nl 

A 2 n/ 2 n 2 

1 — 1 

^ 2 nhn 2 ■ ■ ■ 

A^-V 2 nl 


(3.2) 


J22j 


\n—l 


122 


in—1 


V n—1 . 


0 

Ar‘/i 2 

22 
/32 


(3.3) 


ArV( 


\ n—1 
^3 


A*^ 

A?/ii 

A ^/21 

A?/31 


vn—1 


An-l 

ArVi 2 

ArV 22 


(3.4) 


Ar^ 


/32 


0 

A?/ii 

A ^/21 

A?/31 


\n—1 
^ 2 n 

0 

ArVi 
ArV 22 

ArV 32 


' 2 nJ 


(3.5) 


12 


0 

A?/i 2 

A^/22 

A^/32 


in —1 


(3.6) 
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{Tn)22 — 


0 

1 

0 

A 

0 

A 2 

0 

1 - 1 

1 

A” 

fll 

/l 2 

Al/11 

Ai/12 

^1/11 

Ai/12 

... ArVii 

ArVi2 

A?/i 2 

/21 

/22 

A2/21 

A2/22 

Ai/21 

A^/22 

... ArV2i 

Ar'/22 

A^/22 

/31 

/32 

A3/31 

A3/32 

A3/31 

Ai/32 

... Arv3i 

ArV32 

A3/32 

/ 2 nl 

/ 2 n 2 

A 2 n/ 2 nl 

A 2 n/ 2 n 2 

AL/am 

^ 2 nhn 2 

... A 2 „^/ 2 nl 

A 2 n V 2 n 2 

^ 2 nhn 2 


For the n-fold Darboux transformation, the transformed potentials are 

[4-l=[/o + i[a3, T„], (3.7a) 

(3.7b) 

If we treat the n-fold Darboux transformation as a generalization of the (n — l)-fold, it 
will be multiplications of n one-fold Darboux transformations. It is easy to prove that these 
multiplications have a determinant representation as mentioned above. 

The nth new solutions of the GNLS-MB system after the n-fold Darboux 

transformation will be 

= ^ + 2z(ti)i2, (3.8a) 

, ,,rp X [“2r](r„)ii(T„)i2 + p*iTn)\2{Tn)\2 ” p(71i)ii(T„)ii]| a=-u), (3.8b) 

dctiTjij 

= T7^[^((^n)ll(7;)22 + (7 ;)i2 (T„)2 i) -p*(r„)i 2 (T „)22 +p(r„)n(T„) 2 l]|A=-c.. (3.8c) 

dctylfij 

In order to satisfy the constraints of Darboux transformations, the following conditions must 
be satisfied 


^2k — ^ 


2 fc-l) 


f2k — 



(t = 1, 2. .... n). 


(3.9) 


By calculating, we get det{Tn) = (A — Ai)(A — A 2 )... (A — A 2 n-i)(A — A 2 n), and (ti)i 2 in the 


Eq. (3.8) is the element of the matrix ti at the cross of the first row and the second column, 
where 


1 f {Qn)ll (Qn)l2^ 

W^\\{Qn)2l (Qn)22j 
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(3.10) 










with 



fll 

/12 

^ifn 

^lfl2 

a;/ii 

a;/i2 . 

. Ar‘/i2 

A?/ii 


/ 2 I 

/22 

V/21 

V/22 

y/21 

Ai/22 

. Ar'/22 

a;/2i 

(Qn)ll 

/ 3 I 

/32 

A 3 / 3 I 

A3/32 

Ai/31 

Ai/32 . 

. Ar‘/32 

a;/3i 


f2nl 

/2n2 

^2nf2nl 

^2nf2n2 

AL/2nl 

^2nf2n2 ■ 

. a;-V2„2 

a;„/2„i 


(3.11) 



/ll 

/12 

Al/11 

Al/12 

a;/u 

a;/i 2 . 

• ArTii 

A?/ii 


/21 

/22 

A2/21 

A2/22 

Ai/21 

Ai/22 

• Arv2i 

a;/ 2 i 

{Qn)l 2 

/31 

/32 

A3/31 

A3/32 

Ai/31 

Ai/32 . 

• ArT3i 

a;/ 3 i 


/2nl 

/2n2 

A2n/2nl 

^ 2 nf 2 n 2 

AL/2nl 

Ai„/2n2 • 

• A”-y2nl 

a;„/ 2 „i 


(3.12) 



fll 

fl2 

Al /11 

Al /12 

a;/ii 

a;/i2 . 

. Ar‘/i2 

Ani2 


/21 

/22 

A 2 / 2 I 

A 2/22 

A 2/21 

A 2/22 

. Ar'/22 

An22 

(Qn)21 

/31 

/32 

A 3 / 3 I 

A 3/32 

Ai /31 

Ai /32 . 

. Ar‘/32 

^3^2 


f2nl 

/2n2 

A2n/2nl 

^2nf2n2 

^2n/2nl 

^2nf‘^ri2 

. a;-'/2„2 

^^2nf2n2 


(3.13) 


/ll 

/l2 

Al/11 

Al/12 

a;/u 

Af/.2 . 

• Ar 

‘Til 

Ani2 

/2I 

/22 

A2/2I 

A2/22 

Ai/21 

Ai/22 

. Ar 

^/21 

Ar22 

/3I 

/32 

A3/3I 

A3/32 

Ai/31 

Ai/32 . 

• Ar 

‘T3I 

^3^2 

/2nl 

/2n2 

A2n/2nl 

^2nf2n2 

Ai„/2nl 

^2nf2n2 ■ 

• Ar 

^/2nl 

Ar/2n2 


So far, we discussed the determinant construction of the nth Darboux transformation of 
the GNLS-MB system. Then we will constructed breather solutions and rogue wave solutions 
through these transformations. 


4. BRIGHT AND DARK BREATHER SOLUTIONS OE GNLS-MB SYSTEM 

In this section, we will focus on breather solutions p and r] of the GNLS-MB system, 
which are derived from the periodic seed solutions through the Darboux transformation. Then 
we can get the explicit bright and dark rogue waves of the GNLS-MB system through Taylor 
expansions of the breather solutions. 
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Considering the nonzero background wave, we can take E = p = ifE, and 77 = 1 as 
the initial seeds, where p = az + bt. Then by substituting the seeds into the spectral problem in 
the Eq. (2.1), and by separating variables and superposition principle, the eigenfunction f 2 k-i 
associated with X 2 k-i is given by 


f2k-l — 



Cl'lp 2 k-lliX 2 k-l', t, z) — C2'4^lk-l2{^2k-l''> 

Ci1p2k - 12 { X2k - l \ t , Z ) + C2 ' lp * 2k - ll { Kk - l ^ 


(4.1) 


where 


i^2k-l{X2k-l', t, z) — 


t^ 2 fc-ll(A 2 fc-i; t, z) 
,'ip 2 k- 12 {X 2 k-l] t, z) 


^g(|p+ic(A2fc-i)) 

,^(ci(A2fc-l) + A2fc-1 + 


Here xl) 2 k-iiX 2 k-i] t, z) is the basic solution of the spectral problem in the Eq. (2.1), with 
o, b, d, z, t & M, Cl, 6*2 e C 

a =r(6^ + - Ub'^d^) - 6^ + 2d‘^ + 2/, 

2 


/ 


2 a; — 6 ’ 

c(A2fc-i) =Ci(A2fc-i)2; + C2(A2fc-i)t, 

ci(A2/c-i) =:{(2A2fc-i — b) + t[— 8A2fc_i + 46A2fc_i + (4(i^ — 2b‘^)X2k-i 

I l3i , f 


Qd^b + C] + 


A2fc-1 + 


}c2(A2/c-i); 


C2(A2fc-i) —Xjd'^ + (A2fc-i + -)'^- 


Without lose of generality, in above identities we consider the case when Ci = C 2 = 1, X 2 k-i = 


— | + 7 / 52 fc-i, such that C 2 (A 2 fc-i)^ = d"^ — Im^{X 2 k-i) G M. Under the condition in the Eq. (3.9), 
by substituting eigenfunctions in the Eq. (4.1) into the Eq. (3.8) with the choice of the Eq. 
(3.2) and the Eq. (3.10), we can construct the nth-order solutions E^'^\ pt”! and When 
n = 1 , here we just give an expression of the breather solution E'f^ with the specific parameters 
d=l ,6 = 2 , a;=|,r=|in the following, 

[ 1 ] ^ /3i cos(^i) + (1 - 2^1) cosh(^ 2 ) + 27/3iE:osinh(^2) *(- 19 ^+ 27 ) 

^ /5i cos ( 6 * 1 ) — cosh ( 6 * 2 ) ^ ’ 

where 


(4.2) 


Ko = ^l- 

2 Wo[(— 64 / 5 ^ -|- 8)2; -|- ( 4 / 3 ^ -|- l)t] 


ei = 
02 = 


4/3? + 1 

8piKo{4pf - 19/3j - 3)z 
4/3? + 1 ■ 
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Note that the trajectory of is defined by 


if Kq < 0, and by 


(—64/3^ + 8 ) 2 ; + (4/3^ + l)t — 0 , 


;3l(4;3^ - 1901 - 3)^2 = 0, 


liKl > 0. When /3i = |, we have > 0 then can get the temporal periodic breather solution 
(Ma breather 32 ). Meanwhile the breather solutions 0^^ and 00 can be constructed with 
the same specific parameters of The dynamical evolution of the Ma breather solutions is 
plotted in Fig. 

After a simple analysis, we can know from the eigenfunction that r has an influence on the 
periodic of the breather solutions except Ma breather solutions. Having constructed a bright 
breather for E and dark breathers for p and rj when Kq 0 0, in the next section, our main 
object is to discuss the construction of the rogue wave solutions of the GNLS-MB system when 
Kl goes to zero. 


5. BRIGHT AND DARK ROGUE WAVES OF GNLS-MB SYSTEM 

In this section, firstly, we will construct the first-order bright and dark rogue waves of the 
GNLS-MB system by using the limit method. This kind of solutions only appears in some 
special regions of distance and time and then will be drowned in one fixed non-vanishing plane. 


Under the condition C\ = C 2 = substituting eigenfunctions in the Eq. (4.1) into the Eq. 


(3.8) with Ai = — I + «/3i, by taking the limit /3i d{d > 0), E^^\ and become rational 

When 2 00, t —7- 00 in the 
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solutions E0, p0 and 00 in the form of rogue waves 
expressions of rogue waves, by calculation, we find that the non-vanishing background plane of 
E0, p0 and 00 has nothing to do with r, i.e. 

2d 


\E0\ ^ d, t 


\b — 2ijj\ ’ 


n, 


w ^ 1. 


(5.1) 


Because the solutions are very complicated, we take d = 1, b = 2, oj = \\a order to display 


the results easily. Then the final forms of the rogue waves will be 


— g2i[—(13T-|-3)z-|-t] ^_j_ 

ki 

pin = 2 ie"'l-<'=>-+ 3 ).+<l(i + 


k\ 


n. 


[11 = 1 + 

r I 


16/^5 


(5.2a) 

(5.2b) 

(5.2c) 


where 


ki = 16si2^ — 128s2zt + 200 + 5, 
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k2 = 20 — 14453^2, 

h = -16(17780r2 - 2500r + 101)z^ + 160s4zt - 20t^ - 5, 

= 144siS32;^ — llb2s2S^z^t + ISOss^t^ — 5(86r + l)z + lOt, 
fcs = 16(1580r2 + 740r - 61)z^ - 160s4zt + 20t^ - 5, 


with 


Si = 1940 t — 196r + 29, S 2 = 5t + 1, S 3 = lOr — 1, 54 = 22t — 1. 


(5.3) 


In the Eq. (|5.2a[), from = 0, and = 0, we get that the maximum amplitude 


41] 


of occurs at 2 ; = 0, t = 0 and is equal to 3. The minimum occurs at 2 ; = 0, t = 
and is equal to 0. We also infer that ^ 1 by assuming z ^ 00 , t —>■ 00 , which gives the 

background plane. 


_ _L \/3 


In the Eq. (5.2b), by the same method we get that the height of the background plane is 


2, for —7- 2, when z 00 , t —?• 00 . The maximum amplitude of |pr^| occurs in the form of 

upper ring curve as 

256siz^ - 4096siS2Z^t + 384(4300r^ + lOOr + 91)z^t^ - 5120s2zf^ + 400t^ 

+32(22340t2 + 4940t - 343)z2 - 1280(49r - l)zt + 520^^ - 55 = 0, (5.4) 

the maximum amplitude is equal to \/5- The minimum amplitude of |pr^| occurs in terms of 
four down peaks, 


ddi 


{zi = 
{Z2 = 
{Z3 = 

(Z4 = 


5-75 

72(10r - 1) 

5 + 75 

72(10r - 1) 

5-75 
72(10r- 1) 

5 + 75 

72(10r- 1) 


, h = 
,t2 = - 
■,h = - 

1 ^4 


(590r - 71) - (58r - 37)75 
36(10r- 1)(375 + 2) 
(590r - 71)+ (58r-37)75 




36(10r- 1)(375-2) 
(590r - 71) - (58r - 37)75 
36(10r- 1)(375 + 2) 
(590r - 71)+ (58r-37)75 
36(10r- 1)(375-2) 

Wi 


). 

). 

). 


and is equal to 0. Meanwhile the extreme of the amplitude \pr ^ \ occurs at z = 0, t = 0 and is 
equal to |. 

In the Eq. (5.2c), we conclude that the height of the background plane is 1 because —?• 1 

when z ^ 00 , t —>■ 00 . The rj^r^ has two upper peaks at coordinates as 


{Z2 = 


5 + 75 

72(10t - 1) 


, h = 


(590r - 71) + (58 t - 37)75 
36(10t- 1)(375-2) 


). 
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(Z4 = 


5 + x/5 




(590r - 71) + (58t - 37)\/5 




72(10t - 1) ’36(10t- l)(3y5-2) 
and the height is equal to \/5. The minimum amplitude of r]^r^ occurs at two down peaks with 


coordinates as 


{zi = 
(^3 = 


5 - \/5 
72(10r - 1) 

5 - \/5 


, h = 


, h = 


(590r - 71) - (58r - 37)\/5 
36(10r- l)(3v^ + 2) 
(590r - 71) - (58r - 37)\/5 






72(10r-l)’ 36(10r-l)(3v^ + 2) 

and is equal to — \/5. There is another extreme of the amplitude which occurs at z = 0, t = 0 
and is equal to — y. 

The first-order rogue wave solutions in the Eq. (5.2) are plotted in Fig. [^from which we 
just find three down peaks in the second sub-figure of the Fig. On account of the direction of 
the observation, the two close sub-peaks of the middle down peak are not distinguished clearly 
from the sub-figure. For the same reason, we only catch sight of one down peak in the third 
sub-figure in Fig. 

By the calculation above, we find that r has an impact on the location of extreme except 
when z = 0, but does not change the value of the extreme. We can get a parallelogram by 
connecting the four extreme points (z*, ti) with i = 1, 2 , 3, 4 sequentially. The change of the 
parallelogram reflects the compression and rotation of the rogue wave. The parallelogram is 
plotted in Fig. |^a) with different r. In order to reflect the change of the rogue wave more 
clearly, we also give the upper ring curve of |pr ^ | in Fig. i*-) with the same r in Fig. |^o). The 
area of the parallelogram is equal to 3 e|]^_ 4 | and its function is plotted in Fig. ^c). 

In order to enrich the type of the nth-order rogues, we can modify Ci and C 2 in the Eq. 


(4.1) as following, 


(5.5a) 

(5.5b) 


— g*C2(A2fc-i) Fj=o h'['^2fe-i-(-|+id)].' 

_ g-*C2(A2fe-l) 7)[A2fc-l-(-|-|-*<i)]J 

where Jj G C. Note that \ 2 k -1 = —\ + id is the zero point of C 2 (A 2 fc-i). 

We can get nth-order rogue waves by letting A 2 fc-i ^ — | -|- id in the nth-order breather 
solutions. But it is difficult to obtain higher-order rogue waves from multi-breather solutions 
generally, because the explicit expression of the higher-order breather is very difficult to cal¬ 
culate. In order to solve the problem, we derive the nth-order rogue waves directly from the 


determinant representations of solutions in the Eq. (3.8) by applying Taylor expansions 33 


Theorem 2. For the n-fold Darhoux transformation, substituting the Eq. (5.5) into the eigen 


function in the Eq. (4.1), letting \ 2 k -1 = + id-\- > 0), by adopting Taylor expansions, 
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the determinant representations of the nth-order rogue waves of the GNLS-MB system are given 
as 


eN = E + 2z(t,i)i2, 


(5.6a) 


Pr 


det{T„ 


[—2rj{Trn)l\{Trn)l2 + P* iTrn)l2iTrn)l2 — piTrn)\liTrn)ll]\\=-u 


(5.6b) 


where 


with 


1 ,/rp Ap{{Trn)l\{Trn)22 + {Trn)l2{Trn)2\) ~ P* (^rn) 12 (^rn) 22 + P(^rn)ll (^rn)2l] | 


T = 

± 7*77, 

- 1 / 

(^rn):! 

(^rn) 12^ 

, tj>\ 

1 

1 (Qrn)ll 

(Qrn)l2 j 

|Vb,2n| \ 

^(77n)21 

{Trn) 22 J 

|W^r.2n| 

Y(Qm)21 

(Qrn)22 / 


^ ^01 

^02 

^11 

^12 


”^n-ll 

^n-12 ^ 



'^02 

^01 

~^12 

h\\ 


'^n-12 

^^11 


bbr2n 

'^01 

’h2 

/1?1 

h\2 


”^n-ll 

”^n-l2 
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^2n—1* 
\ '^02 

^2n—1* 
'^01 

^2n—1* 

'n2 

L2n—1* 

’hi 

^2n—1* 
'^n-12 

^2n-U , 

’h-ll / 



1 

0 

A 

0 


An-l 

0 

A” 


^01 

^02 

^11 

^12 


^n-11 

Ki-12 

^nl 

{Trn)ll = 

~^02 

'^01 

"12 

Ml 


"n-12 

hU 

’’^n-ll 

-^i*2 

hh 

’h2 

/I?! 

/^?2 


”^71-11 

”^71-12 

^nl 


L2n—1* 
'^02 

^2n—1* 
'^01 

L2n—1* 

'n2 

L 2n—1* 


^2n—1* 
'^n-12 

L2n—1* 

L2n—1* 

'V2 


0 

1 

0 

A 


0 

An-l 

0 


’hi 

^02 

^11 

^12 


^n-ll 

^n-12 

’hi 

{Trn) 12 = 

'^02 

"01 

-^2 

Ll* 

"11 


"n-12 

hi* 

"n-11 

^7i2 

’hi 

’h2 

’hi 

’h2 


”^n-ll 

'^n-12 

’hi 


^2n—1* 
'^02 

^2n—1* 
'^01 

^2n—1* 

'n2 

^2n—1* 


^2n—1* 
'^n-12 

^2n—1* 

^2n—1* 

’h2 


X=—(jj 5 

(5.6c) 


(5.7) 
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(^Qrn)22 


1,1 

'^01 

'^02 

Ml 

M 2 



hl2 

~b^2 


”^12 

h\l 


~^ n-12 


hh 

r3 

'^02 

Ml 

M 2 


r3 

hl2 

7 2n — 1 * 
'^02 

T2n — 1 * 
^01 

7 2n — 1 * 
^12 

T2n — 1 * 
^11 


T2n — 1 * 
'^n-l2 

7 2n — 1 * 

^nl 

ft' = A 

mj Q^l 

[i-^+id+E^r 

MM = 

b 

~2 

id 


m = 

0 , 1 , 2 , • 

• • , n; 

3 = 1 , 2 ; 

1 

= 1 , 2 , •••, 

2n. 


According to our analysis, we can know that there are n + 4 free parameters denoted as 
(Jo, Ji, • • • , Jn-i, b, d, u, r) in the nth-order rogue wave solutions. Next, we will consider the 
types of the nth-order rogue waves with these parameters. For convenience, we let d = 1, 6 = 
0 , a; = |, r = I in the following. 


For n = 2, the Eq. (5.6) in Theorem 2 can give the second-order rogue wave solutions of 


the GNLS-MB system, when Jq = Ji = 0, we can get the fundamental pattern which is plotted 
in Fig. 1^ the maximum amplitude of \E^^\ is five times as high as the background plane 


34 


When Jo = 0, Ji = 100, the fundamental pattern can be split into three first-order rogue waves 
(triangular structure [^), which is shown in Fig. 


For n = 3, the Eq. (5.6) can give the third-order rogue wave solutions of the GNLS-MB 


system. When Jo = 0, Ji = 0, J 2 = 0, we can get the fundamental pattern as plotted in Fig. 
1^ the maximum amplitude of is seven times the height of the background plane. When 
Jo = 0, Ji = 100, J2 = 0, the fundamental pattern is split into six first-order rogue waves 
(triangular structure), which is shown in Fig. But when Jo = 0, Ji = 0, J 2 = 5000, we get 
the ring structure which is also made up of six first-order rogue waves as in Fig. Note that, 
the ring structure and the triangular structure have the same profile when n = 2. 

So far, we have obtained three types of structures: fundamental pattern, triangular struc¬ 
ture and ring structure. According to the conclusions and decomposition rule of the higher-order 


rogue waves in Ref. |M,36 , we know that higher-order rogue waves are generally the combina¬ 


tion of the above three types. The specific application of these conclusions will be given in the 
following. The parameters d, b, uj, t take the same values as above. 


For the nth-order rogue waves in the Eq. (5.6), it will result in the fundamental pattern if 


the parameters Jj are all zero, then the fundamental pattern can split into triangular structure 
when Jj = 0 except Ji. The parameter J^_i determines the ring structure when n > 3, 
Jn-i ^ 1 and other parameters Jj are all zero. The ring structure consists of a fundamental 
pattern of (n — 2)th-order rogue waves located in the center and 2n — 1 first-order rogue waves 
located on the outer circle. As Eig. |^when n = 4 with parameters J = 1, 6 = 0, a; = |, r = 
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Jo = 0, Ji = 0, J2 = 0, J3 = 10®, their corresponding density plots are shown in Fig. 

The fundamental pattern of {n — 2)th-order rogue waves in the center can keep on splitting into 
triangular structure or ring structure by choosing the parameters Jj(j = 0, 1, • • • , n — 3). And 
the process can go on indehnitely until they completely decompose. For example, we can get 
the ring structure of the fifth-order rogue waves by setting J4 7^ 0. In this case: if we assume 
Ji 7^ 0, the inner third-order rogue waves can split into triangular structure (Fig. 0 ; if we 
set J2 7^ 0, the inner can split into a ring structure (Fig. [I^ . Other parameters in addition to 
Jn-i can also generate rings which are different from the above. For example, if and only if the 
parameter J„_2 7^ 0, the nth-order rogue waves have a double-ring structure, which has two 
rings with 2n — 3 first-order rogue waves for each ring and a fundamental pattern of (n — 4)th- 
order rogue waves (Fig. [^. And also if and only if J„_3 7^ 0, we can get a decomposition 
with three rings (Fig. [^. We can also create more composite structures by selecting different 
combination of parameters appropriately. Here we will no longer show them in detail. 

In fact, Jo has no influence on the type of the nth-order rogue waves. It can just adjust 
the position of the nth-order rogue waves on the background plane 37 . Next, we will discuss 
the nonlinear superposition of the rogue wave and breather solutions. Under the condition 


Cl = C 2 = 1, X 2 k-i = ~| + iP 2 k-i, substituting eigenfunctions in the Eq. (4.1) into the Eq. 
(3.8), we can get the nth-order breather solutions. We know that the nth-order rogue waves 
could be generated by assuming P 2 k-i ^ d{k = 1, 2, • • • , n). When n > 2, if we only take the 
limit for I52k-i{k = 1, 2, • • • , n^), 1 < n^ < n, we can construct new solutions which are the 
nonlinear superposition of n^th-order rogue waves and (n — n^jth-order breather solutions. In 
order to separate the n^th-order rogue waves from the (n—n^jth-order breather solutions, we can 
modify Ci and C 2 by substituting the Eq. (5.5) into eigenfunctions in the Eq. (4.1) only when 
A: = 1, 2, • • • , Note that, we take (Ji = (72 = 1 in the Eq. (4.1) when k — rir + 1, ■ ■ ■ , n. 

Due to the tediousness of the solutions, we give the dynamical evolution of the results. 
For n = 2, Ur = 1, when the two components are separated completely, the solutions appear 


as a first-order rogue wave and a first-order breather solution as plotted in Fig. 15, which are 


similar to Fig. 3(a) in Ref. 38 . However, when the two components are nonlinearly overlaid, 
the solutions appear as a first-order breather with a central fundamental pattern of a second- 


order rogue wave as plotted in Fig. 16, which are similar to Fig. 3(b) in Ref. 38 . Comparing 
Fig. [T^ with Fig. we observe that the three first-order peaks merge into a second-order 
peak, in other words, a second-order peak splits into three first-order peaks. These hybrid 
solutions give the interaction between the rogue wave and breather solutions. 
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6. CONCLUSION 


In this article, we derived the Darboux transformation of the GNLS-MB system which 
describes the ultrashort pulse propagation in the resonant erbium-doped nonlinear optical fiber 
with higher-order effects, and constructed the determinant representation of the n-fold Darboux 
transformation in Theorem 1. We established the determinant representation of the solutions 
(£;W, generated from the known trivial seed solutions {E, p, rj) in the Eq. (|3.8|). 

Under the reduction conditions \ 2 k = ^ 2 k-i hk = ( the breather solutions 

from the nonzero periodic seed solutions. Using the method of limit technique and Taylor 
expansions, we constructed the determinant representation of the nth-order rogue waves in 
Theorem 2. The rogue waves show interesting characteristics that might attract physicists in 
experiments to observe them, in contrast with the common bright rogue wave the dark rogue 
waves for p and 77 have two (or more) dominant down peaks in their profiles, and there is an 
upper ring in the profile of p. 

By analyzing higher-order terms, we found that r has an effect on the periodic of the 
breather solutions except Ma breather solutions. By calculating, we knew p and 77 had the 
same locations in connection with r and that r did not affect the values of their extremes. 

In the last section of this article, we also gave the three basic types of the Trth-order rogue 
waves, and combination structures could be obtained by choosing proper parameters. At the 
end of the last section, we introduce the hybrid solutions which are the nonlinear superposition 
of the rogue wave and breather solutions. These solutions display the interaction between 
the rogue wave and breather solutions, which could help us to understand the generation of 
rogue waves better. Moreover, the exploration for other coupled systems with more complex 
higher-order optical effects can be done in our future work. 
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Figure 1. (Color online) The first-order breather solutions (F, p, 77) of the 
GNLS-MB system when d = 1, 6 = 2, a; = |, r = |, /3i = |. 



Figure 2. (Color online) The first-order rogue wave solutions {E, p, rf) of the 
GNLS-MB system when d= 1, b = 2, u) = r = 




Figure 3. (Golor online) When d — 1, b = 2, co = ^: (a) the parallelogram 
consists of the four extreme points (zi, ti) with r = 0 (blue, solid) and t = \ 
(red, dash). (6) the upper ring curve of the |pr^| with r = 0 (blue, solid) and 
r = I (red, dash), (c) the area function of the parallelogram in (a) with different 

T. 
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Figure 4. (Color online) The fundamental pattern of the second-order rogue 
wave solutions {E^p^rj) of the GNLS-MB system when (i=l,6 = 0,a;=|,r = 
Jo = O5 = 0- 



Figure 5. (Color online) The triangular structure of the second-order rogue 
wave solutions {E^p^rf) of the GNLS-MB system when d = l ,6 = 0,cj = |,'r = 
Jo = 0, Ji = 100. 



Figure 6. (Color online) The fundamental pattern of the third-order rogue wave 
solutions (£”, p, 77 ) of the GNLS-MB system when (i=l,6 = 0, cj=|,r = 
Jo = 0 , Ji = 0 , J2 = 0 . 
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Figure 7. (Color online) The triangular structure of the third-order rogue wave 
solutions (F, p, 77) of the GNLS-MB system when d = 1, b = 0, u = t = 
Jo = 0, Ji = 100, J 2 = 0. 



Figure 8 . (Color online) The ring structure of the third-order rogue wave solu¬ 
tions {E, p, r]) of the GNLS-MB system when d = 1, b = 0, oj = ^, r = ^, Jq = 
0, Ji = 0 , J 2 = 5000. 



Figure 9. (Color online) The ring structure of the fourth-order rogue wave 
solutions {E, p, 77) of the GNLS-MB system when d = 1, b = 0, lo = t = 
Jo = 0, Ji = 0, J 2 = 0, J 3 = 10®. 
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Figure 10. (Color online) Density plots of the wave amplitudes for 
ri^r^) when d = 1, b = 0, cu = t = Jo = 0, Ji = 0, J2 = 
0, J3 = 10®. They have a fundamental pattern in a ring. 



-3-2-10 1 2 3 
z 




z 


Figure 11. (Color online) Density plots of the wave amplitudes for 
(iFr^'l, rir^) when J = 1, 6 = 0, n; = |, r = |, Jq = 0, Ji = 100, J2 = 

0, J3 = 0, J4 = 10®. They have a triangular structure in a ring. 
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Figure 12. (Color online) Density plots of the wave amplitudes for 
(IF^I, when d = 1, b = 0,n; = ^, r = |, Jo = 0, Ji = 0, J2 = 

5000, J3 = 0, J4 = 10®. They have a ring structure in the outer ring, the outer 
ring is made up of nine first-order rogue, and the inner ring has hve first-order 
rogue waves. 
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Figure 13. (Color online) Density plots of the wave amplitudes of 
(iFi-^'l, 7]^) when J = 1, 6 = 0, a; = ^, r = |, Jo = 0, Ji = 0, J2 = 

0, J3 = 10®, J4 = 0. They have two concentric rings, each of them consists of 
seven first-order rogue waves. 
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Figure 14. (Color online) Density plots of the wave amplitudes of 
(IF^I, 7^^) when d = 1, 6 = 0, a; = r = Jo = 0, Ji = 0, J2 = 

5000, J3 = 0, J4 = 0. They have three concentric rings, each of them consists of 
five first-order rogue waves. 



Figure 15. (Color online) The completely isolated form of the solutions which is 
the nonlinear superposition of the first-order rogue waves and first-order breather 
solutions, when d=l, 6 = 0 , c<;=|,r = |,/33 = |, Jo = — lOh 



Figure 16. (Color online) The form shown as a first-order breather with a 
central fundamental pattern of second-order rogue waves, when d=l,6 = 0,c<; = 
I, T =1, P3= |, Jo = 0. 
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